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SUMMARY 



This paper is the first of a series covering an in- 
vestigation of a family of airfoils all formed from a 
"basic profile. It gives in preliminary form the results 
obtained from tests in the N.A.C.A. variable density wind 
tunnel of six symmetrical airfoils, differing only in max- 
imum thickness. The maximum, thickness-tu-chord ratios are 
0.06, 0.09, 0.12, 0.15, 0.18, and 0.21. The results are 
analyzed with a view to indicating the variat ion of the 
aerodynamic characteristics with profile thickness. 



INTRODUCTION 



The. forme of the airfoil sections that are in common 
use ; today are the result of a more or less systematic in- 
vestigation made at G-ottingen of a large number of air- 
foils. Previously airfoils such as the. R.A.E. 15 and 
U.S.A.' 27, developed from airfoil profiles investigated 
in England, were widely used. Because most airfoils have 
been developed from low-scale tests, the forms developed 
may not be the optimum for full-scale values of the 
Reynolds Number. A number of airfoils have been investi- 
gated in the variable density wind tunnel at values of 
the Reynolds Number approaching, those of flight (refer- 
ence l), but with the exception of the M series and a 
series of propeller sections, the airfoils have not been 
related in such a way that the results could be satis- 
factorily correlated. 

The object of an investigation now being carried out 
by the National Advisory Committee for Aeronautics is to 
obtain the characteristics, at large values of the Reynolds 
Number , of a wide variety of related airfoils. The bene- 
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fits of a systematic investigation of airfoil profiles at 
large values of the Reynolds lTunber are so self-evident 
that it is ..hardly necessary to point then out. Hot only 
do the results of such investigations greatly facilitate 
the choice of the most satisfactory airfoil for a given 
application but, "because the results may be correlated to 
indicate the trends of the aerodynamic characteristics 
with changes of shape, they may point the way to the de- 
sign of new shapes having better characteristics. 

Airfoil profiles may be" considered as made up of 
certain profile thickness forms disposed about certain 
mean camber lines. The major shape variables then be- 
come two: the thickness form and the mean camber line 
form. The thickness is of particular importance from a 
structural standpoint. On the other hand, the form of the 
mean camber line determines almost independently some of 
the most important aerodynamic properties of the airfoil 
section, e.g. , the pitching moment characteristics and the 
angle of zero lift'. 

The related airfoil profiles, were obtained for this 
investigation by changing systematically these shape vari- 
ables. A single basic thickness variation was chosen for 
the first group of airfoils. Sections having a different 
maximum thickness were obtained by the application of a 
factor to all the basic ordinates. The following ratios 
of maximum thickness-to-chord were chosen: 0.06, 0.09, 
0.12, 0.15, 0.18, and 0.21. The cambered profiles were 
then obtained by combining these thickness forms with dif- 
ferent mean' camber lines. Since this report does not deal 
with the cambered airfoils, it will be sufficient to state 
that the various mean camber line forms are obtained by 
varying the maximum camber and by varying the distance 
from the loading edge to> the position of the maximum cam- 
ber. The airfoils so produced are designated by a numbor 
of four digits; the f-irst indicates the maximum mean cam- 
ber; the se,cond, the position of the maximum mean camber; 
and the last two, the maximum thickness. Thus the N.A.C.A. 
2315 airfoil has a maximum mean camber of 2 por cent of 
the chord at a position 0.3 of the chord from the leading 
edge, and a maximum thickness of 15 per cent of the chord; 
the U.A.C.A. 0012 is a symmetrical airfoil having a maxi- 
mum thickness of 12 per cent of the chord. 

This note presents in preliminary form the results of 
the tests of the first group of six airfoils. These tests 
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were mado in April, 1931. Since these airfoils may "be con- 
sidered as the "basic soctions from which the others will he 
formed, the results are of particular importance. It is 
therefore considered desirable to present these results "De- 
fore the tests of all the airfoils of the series are com- 
pleted and the results analyzed. 



DESCRIPTION 02 AIRFOILS 



Well-known airfoils of a certain class, including the 
Gottingen 3S8 and the Clark Y , which have proved to he ef- 
ficient, are nearly alike when their mean camhor is removod 
and they are reduced to the same thickness. A thickness 
variation similar to that of these airfoils was therefore 
chosen for the development of the S.A.O.A. airfoils. A 
formula defining the shape was used as a method of produc- 
ing fair profiles. 

If the chord is taken along the x axis from 0 to 1, 
the ordinates y are given "by a formula of the type 

±y = a 0 & + a x x + a 3 x 3 + a^x 3 + a 4 x 4 . 

The equation was adjusted to" give the desired shape "by im- 
posing the following conditions to determine the constants: 

^ 1. Maximum ordinate 0.1 at 0.3 chord **" 

y = 0.1 at x = 0.3 

? = 0 at I = 0.3 
dx 

2. Thickness of trailing edge <■ 

y = 0.002 at x = 1 

3. Trailing edge angle 

^ = -0.234 at x = 1 
dx 

4. Nose shape 

y = 0.078 at x = 0.1 
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The ordinates of the "basic section, which is a symmetrical 
section having a total thicknoss of twenty per cent of the 
chord., are then given "by the formula 

±7 = 0.296900 Jx - 0.126000X - 0.351600x s + 0.284300x 3 
- 0.101500 x * 

The curve of the "basic section, shown in Figure 1, 
may "be compared with the plotted points on the same figure 
that were obtained "by removing the mean cam"ber from the 
G-Bttingen 398 and the Clark Y sections, and applying a 
factor to the resulting thickness curves to bring them to 
the same maximum thickness. 

The ordinates of the basic section as given by the 
formula are multiplied- by a factor to obtain the related 
sections having any desired maximum thicknoss. The lead- 
ing-edge radius for the basic section is found from the 
formula to be 




or, for a derived section of maximum thickness-to-chord 
ratio t, the leading-edge radius expressed as a ratio 
to the chord is given by 

= (0.2969 ) 2 t a f 
0.8 

= 1.10 t 3 

The ordinates of the six besic sections having thicknoss- 
to-chord ratios 0.06, 0.09, 0.12, 0.15, 0.18, and 0.21 are 
given in Table I. 

The airfoil models are made of heat-treated duralunin. 
A special airfoil generating machine is employed that works 
from a six-fold templet of the section. The templets are 
carefully laid out on a layout table that permits the plot- 
ting of the stations and ordinates to an accuracy of 0.001 
inch. The templets are then cut out and checked for pre- 
cision of contour. A section of the model airfoil is also 
checked after the cut is started and any necessary correc- 
tions are made on the templet. As the cut progresses, the 
maximum thickness is checked from time to tine to guard 
against errors resulting from excessive tool woar . 
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The models are hand-finished .to remove small, tool 
marks and then "buffed to produce a polished surface. To 
insure accuracy of alignment in the tunnel, a special 
drilling jig is employed to drill the airfoil models for . 
mounting. 



TESTS 



The tests were conducted -in the redesigned variable 
density, tunnel . A report describing the tunnel and the 
details' of airfoil testing is "being prepared "but, "because 
this report is not yet available and "because the tunnel 
has "been m'odified in several important respects since a 
description of it has "been published, a short description 
of the tunnel and of the details of tho tests will be in-, 
eluded here. Several changes will be noted if the dia- 
grammatic section of the tunnel shown in Figure 2 is com- 
pared with previously published ones. (See references 2 
and 3.) As a result of these changes the air flow in the 
tunnel has been improved. The variation of the air-flow 
direction, as indicated by a yaw head when passed across 
the test section, is less than ± l/4°. The velocity 
distribution at the test section is also very satisfac- 
tory. Aside from a small central region, approximately 
3 inches in diameter, in which the velocity is less than 
1 per cent low, the dynamic pressure is approximately 
uniform across the test section. 

The models are mounted in the tunnel .by. means of pins 
at the upper ends of the two main supporting struts shown 
in Figure 2. The pins are located, with reference to the 
model, on the chord line one-quartor of the chord behind 
the leading edge; and with reference to tho balance, in 
line with the intersection of the drag and lif^ balance 
linkages. (See fig. 2.) The location of the pins with 
respect to the balance is checked from time to time and 
changed, if necessary, so that loads applied to the pins 
produce no deflection of the 'moment balance. The accuracy 
of the moment balance has also been increased, so that ac- 
curate moment values measured about the quarter chord 
point may be obtained directly from the balance. 

The airfoil tests wer6 made at only one : value of the 
Reynolds llumbor, approximately 3,000,000. This value was 
obtained by using an air pressure in the tunnel of approx- 
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imately 20 atmo spheres •( the highest at which the tunnel 
is usually ' operated) ,■ and "by using_a slightly reduced 
speed to maintain the ..dynamic pressure at such a value 
that the c ount ers ,. which read the balance leads, indicated 
directly the values of the coefficients. It was thus pos- 
sible to plot the final curves of the coefficients as the 
test progressed. The Reynolds Number of 3,000,000 corre- 
sponds approximately to. that reached by most airplanes in 
flight near minimum speed. The method of testing was 
otherwise nearly the same as that described in reference 3. 

'The tare drag was determined by measuring the drag of 
the supporting members while they were - connected inside a 
hollow dummy airfoil mounted independently of the balance. 
These measurements were made with the dummy airfoil at sev- 
eral angles of- attack. The tare drag at zero angle, ex- 
pressed as a coefficient based on the wing area, was found 
to be 0.0055. 

The zero angle attitude of the airfoil and the drag 
coefficients were corrected for air flow misalignment. A 
slight upflow in the neighborhood of the airfoil results 
from, the damming effect of the support struts and support- 
strut fairings. The effective value of the upflow angle, 
as determined from tests of airfoils in erect and inverted 
positions, is 0.0050 radian. One-half of one per cent of 
the lift is therefore added to the drag and the zero angle 
is set slightly below the horizontal. The alignment of 
the balance with the horizontal is checked from time to 
time by noting the drag balance deflection as weights are 
placed on the main balance frame. 



RESULTS 



The results are presented in tabular and in graphic 
form. The lift coefficient, angle of attack, drag coeffi- 
cient and moment coefficient are given in Tables II to VII. 
The lift coefficient C L is considered as the independent 

variable. The lift coefficient values as measured in the 
tunnel are therefore given in the tables. The angle of 
attack a 0 is the angle of attack' corrected to infinite 

aspect ratio by the method described in reference 1. The 
angle of attack for aspect ratio 6. may be obtained by 
adding 3.58 degrees. The drag coefficient C D is the 
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profile drag coefficient obtained from the test data by 
deducting the induced drag coefficient calculated by the 
method described in reference 1. The drag coefficient for 
aspect ratio 6 may he found by adding 0.0559 C^ 2 . The 

moment coefficient Cm , ie the coefficient of the pitch- 

™C/4 

ing moment a'bout the quarter chord point and may he con- 
sidered as independent of aspect ratio. 

The profile drag coefficients for all the sections 
are plotted against the lift coefficient in Figure S. The 
lift coefficient curves for all the airfoils are plotted 
in Figure 4. 



DISCUSSION 

Variation of the Aerodynamic Characteristics 

with Thickness 



The minimum drag, coefficients for the airfoils are 
plotted against thickness in Figure 5. It will he noted 
that the minimum drag coefficient increases progressively 
with thickness. The minimum profile drag coefficient 
C Dq minimum may he represented as a function of the thick- 
ness by 

Cr> \ min, = 0,0065 + 0.0083t + 0.097t 2 

where t is the maximum thickness expressed as a fraction 
of the chord. The above function is represented as the 
curve in Figure 5. The points on. the same figure are min- 
imum profile drag values taken, from the faired profile 
drag curves* 

It is of interest to consider the minimum drag coef- 
ficients based on projected frontal area. This type of 
coefficient, which is employed when the sections are used 
as strut sections, nay.be designated C D '. By dividing the 

above expression by t, 

C D ' = 0^0^5 + 0.0083 + 0.0972t 
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_D_ 0^00 6.5 
dt t a 



+ 0.0972 



and equating to zero, 



t 



a 0.258 



Cy,' min. = .0586 



In other words, the coefficient "based on frontal area de- 
creases with increasing thickness throughout the thickness 
range investigated, and would apparently continue to de- 
crease until a thickness of about 26 per cent was reached. 

The variation of the maximum lift coefficient with 
thickness is indicated in the following table: 



It will "be noted that the highest maximum lift coefficients 
are obtained with the moderately thick airfoils. It should 
bo noted, however, that the lift curves (fig. 4) which have 
very high maximum values have also comparatively sharp 
peaks so that in practice the highest maximum lift coeffi- 
cients may be of little or no value. Subsequent tests of 
the 0012 airfoil have indicated that the flow is so critic- 
al in the region of maximum lift that the results are dif- 
ficult to reproduce. The estimated uncertainty of the val- 
ues given for the maximum lifts may be as large as 5 per 
cent for the airfoils having lift curves that indicate an 
abrupt break boyond the maximum. 

The lift-curve slope changea with thickness as shown 
in Figure 6. The points on the figure represent the de- 
duced slopes for an infinite span wing. It will be noted 
that all of the values lie below the theoretical value 
for thin wings, 2W ; per radian. These results are in ac- 
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cord with previous results (reference l) in that the' lift- 
curve slope tonds to decrease with increasing thickness. 

The moment coefficient for zero angle of attack would 
be zero for all the airfoils if the flow and the airfoils 
were exactly symmetrical. Thin-airfoil theory indicates 
that the pitching moment about a point one- quarter of the 
chord behind the leading edge is zero. Actually, these 
results indicate that the. pitching moment . is ;very nearly 
zero at zero angle, but that the moment about the quarter 
chord point, instead of remaining zero, increases with the 
angle of attack, as shown in Figure 7, until the angle. of 
attack approaches that of maximum lift. The .'departure 
from 5 the the-ory increases with increasing thickness. This 
subject will be discussed further under.. the heading, 
"Variation of Characteristics with Lift or Angle." 

The general efficiency of an airfoil can not be in- 
dicated by means of a single number.' ' The ratio of the max- 
imum lift to minimum profile drag is, however, of some val- 
ue as a measure of tho efficiency of an* airfoil ' section. 
This ratio is largest for the sections between. 9 and 12. per 
cent thick and falls off rapidly when the thickness is in- 
creased beyond 18 per cent. The variation of this ratio 
with thickness is indicated. in Figure 8. 

Variation of Characteristics with Lift or Angle 



The variation of tho profile drag coefficient with 
lift coefficient for all the sections is shown in Figure 3. 
In reference 4 the variation of the profile drag coeffi- 
cient with lift coefficient was approximately represented 
by a single function of the lift coefficient for all air- 
foils. It is evident, however, that the present results, 
which are more accurate and cover a greater range of thick- 
ness, do not indicate the same variation of the drag with 
lift for the different airfoils. 

If the increase of the profile drag coefficient with 
lift coefficient could be approximately represented as 
being proportional to the square of the lift coefficient, 
then, following the method developed for performance, pre- 
diction at California Institute of Technology, the addi- 
tional profilo drag could be included with the induced 
drag. The additional profile drag coefficients 
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Ctj - Ct) min. are plotted against the square of the lift 
coefficient in Figure 9. The straight line 

C D 0 " C D Q min * " 0.0062C L 2 

represents the additional profile drag coefficient of the,. 
0012, 0015., and 0018 airfoils for any lift coefficient be- 
tween 0 and 1 to a precision of ± 6 per cent of the mini- 
mum drag of the 0015. Combining the above equation with 
the empirical equation previously • developed to represent 
the variation of the minimum profile drag coefficient with 
thickness, the following equation is obtained. Lt repre- 
sents the profile drag coefficient of the moderately thick 
N.A.C.A. symmetrical airfoils at .values of the lift coef- 
ficient below 1 . 

°D 0 = 0.00S5 + 0.0083t + 0.0972t 2 + 0.0062C L S 

It is of interest to note that the factor 0,0Q62C L 2 is 

11.7 per cent of the induced drag of an elliptical wing 
of aspect ratio 6. 

The fact that the moment coefficients near zero lift 
increase with increasing lift coefficients (fig. 7) indi- 
cates that the center of pressure is' ahead of. the 25 per 
cent chord point. The following table indicates the 
points about which the moment coefficients are constant. 
In other words, the airfoils are stable in pitch for mod- 
erate angles of pitch about axes ahead of this point and 
unstable about axes behind the point. 



Airfoil Distance from l/4 chord point 

forward to point of zero moment 
(per cent of chord) 

1T.A.C.A. 0006 0.8 

N.A.C.A. 0009 .7 

1I.A.C.A. 0012 1.2 

H.A.O.A. 0015 1.4 

N.A.C.A. 0018 1.7 

N.A.C.A. 0021 2.1 
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The centers of pressure for small angles of attack 
are shown to "be farther forward for the thick airfoils. 
The differences, however, are small and are, so influ- 
enced by tip effects that the equilibrium points indi- 
cated in the table may not be taken as accurately rep- 
resenting properties of the sections. 



These results show that the aerodynamic character- 
istics of related symmetrical airfoils vary systematical- 
ly with the thickness ratio of the section. The highest 
value of the ratio of maximum lift to minimum drag was' 
found to. correspond with a thickness ratio of a little 
less than 0.12. 
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Ordinates of N.A.C.A. Airfoils 



Station 






Ordinates 


chord) 


% chord 


has ic 


0006 


uuuy 


OOi<2 


OOlO 


0018 


AAA1 

0021 


n 




n 


0 


0 


0 


0 


0 


i ? 

X « Ci u 


"7 T c e C 
O • XuOO 


Q A"7 


1.420 


1.894 


2.367 


2.841 


3.314 




rb . o O f 5? 


t mi 

-L . OU f 


1.961 


2.615 


3.268 


3.922 


4.576 






X . f f t 


2.666 


3.555 


4.443 


5.332 


6.221 


7.5 


6.9998 


2.100 


3.150 


4.200 


5.250 


6.300 


7.350 


10 0 




9 *X AT 

d « CtI 


3.512 


4. 683 


5.853 


7.024 


8.195 


15 


8.9086 


2.673 


4.009 


5.345 


6.681 


8.018 


9.354 


20 


9.5626 


2.869 


4.303 


5.738 


7.172 


8.606 


10.041 


30 


10.0029 


3.001 


4.501 


6.002 


7.502 


9.003 


10.503 


40 


9. 6717 


2.902 


4.352 


5.803 


7.254 


8.705 


10.155 


50 


8.8234= 


2.647 


3.971 


5.294 


6.618 


7.941 


9.265 


60 


7.6056 


2.282 


3.423 


4.563 


5.704 


6.845 


7.986 


70 


6.1065 


1.832 


2.748 


3.664 


4.580 


5.496 


6.412 


30 


4.3719 


1.312 


1.9 67 


2.623 


3.279 


3.935 


4.590 


90 


2.4128 


.724 


1.086 


1.448 


1.810 


2.172 


2.533 


95 


1 . 3443 


.403 


.605 


.807 


1.008 


1.210 


1.412 


100 


.2100 


.063 


.095 


.126 


.158 


.189 


.221 


I .E.Rad. 


4.405 


.394 


.887 


1.576 


2.464 
,. . 


3.549 


4.830 
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SABLE II 

Airfoil: N.A.G.A. 0006 

Average Reynolds Number: 3,120,000. 
Size of model: 5 X 30 inches. 
Pressure, standard atmospheres: 20.8. 
Test No.: 555. Variable Density Tunnel. 
Date: April 8, 1931. 









C M / 

C/ 4. 


-0.847 


-11.3 


0.1582 


0.055 


-.786 


-9.5 


.1324 


.020 


-.613 


-6.0 


.0129 


-.005 


-.311 


-3.0 


.0078 


-.003 


.004 


0 


.0070 


0 


.157 


1.5 


.0077 


.002 


.315 


3.0 


.0084 


.003 


.620 


6.0 


.0135 


.004 


.794 


9.5 


.1195 


-.025 


.848 


11.3 


.1645 


-.059 


.867 


13.2 


.2125 


- .090 


.857 


15.3 


.2528 


-.110 


. 8.35 


17.3 


.2877 


-.125 


. 823 


21.4 


.3503 


- .135 


.819 


27.4 


.4695 


-.145 
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TABLE III 

Airfoil: H.A.C.A. 0009 

Average Reynolds Number: 3,110,000. 
Size of model: 5 X 30 inches. 
Pressure, standard atmospheres: 20.6 
Test Ho.: 558. Variable Density Tunnel. 
Date: April 9, 1931. 





a 0 




Cm c/ 4 


-1.048 


-10.7 


0.0175 


-0.005 


-.905 


-9.1 


.0137 


-.003 


-.602 


-6.1 


.0101 


-.004 


-.298 


-3.1 


.0083 


-.002 


.011 


0 


.0080 


.001 


.165 


1.5 


.0082 


.003 


.316 


3.0 


.0088 


.001 


. 622 


6.0 


.0104 


.004 


.915 


9.1 • 


.0144 


.003 


1.064 


10. 6 


.0170 


0 


1.195 


12.2 


.0312 


-.001 


1.099 


14.5 


.1383 


-.049 


1.037 


13.7 


.2255 


-.085 


.909 


21.1 


.3417 


-.121 


.813 


27.4 


.4614 


-.138 
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TABLE 17 

Airfoil: H.A.C.A. 0012 

Average Reynolds Humher: 3,230,000. 
Size of model: 5 X 30 inches. 
Pressure, standard atmospheres: 20.4. 
Test Ho.: 562. YariaDle Density Tunnel. 
Date: Ax>ril 13, 1931. 



°L 


a 




C M / 








C/4 


-1.187 


-12.2 


0.0194 


-0.002 


-.904 


-9.1 


.Q137 


-.003 


-.604 


-6.1 


.0108 


-.005 


-.289 


-3.0 


.0090 


-.002 


.003 


0 


.0089 


-.001 


.155 


1.5 


.0090 


.001 


.310 


3.0 


.0096 


.002 


.610 


6.1 


.0111 


.008 


.912 


9.1 


.0141 


.011 


1.195 


12.2 


.0207 


.008 


1.325 


13.8 


.0252 


.003 


1.413 


15.2 


.0337 


.003 


1.160 


16.9 


.1562 


-.054 


.990 


20.9 


.2811 


-.096 


.864 


27.3 


.4295 


-.128 



t 
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TABLE V 

Airfoil: 1T.A.C.A. 0015 

Average Reynolds ITumber: 3,110,000. 
Size of model: 5 X 30 inches. 
Pressure, standard atmospheres: 20.4. 
Test Ho.: 530. Variable Density Tunnel. 
Date: April 10, 1931. 





a 

0 


°»o 
o 


C/ 4 


-1.195 


-12.2 


0.0204 


-0.002 


-.915 


-9.1 


.0145 


-.001 


-.611 


-S.l 


.0119 




-.312 


-3.0 


.0102 


-.001 


-.009 


0 
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TABLE VI 

Airfoil: tf.A.C.A. 0018 

Average Reynolds Number: 3,140,000. 
Size of model: 5 X 30 inches. 
Pressure, standard atmospheres: 20.8 
Test Ho.: 553. Variable Density Tunnel. 
Date: April 7, 1931. 
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TABLE VII 

Airfoil: N.A.C.A. 0021 

Average Reynolds Number: 3,180,000. 
Size of model: 5 X 30 inches. 
Pressure, standard atmospheres: 20.6. 
Test Ho.: 559. Variable Density Tunnel. 
Date: April 10, 1931. 
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ffig.l Thickness variation. 




Pl g .2 Diagranmtle section of tae reaesigxww Taruwj. 
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ITig.H Lift curves corrected to infinite aspect ratio. 
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Fig.!3 Variation of ndninurn profile drag with, thickness. 
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2 1*. 6 S 10 12 lH 16 IS 20 22 
Maximum thickness , $ chord 

Fi^.6 Variation of lift-curve slope with thickness. 
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Pig.7 Moment coefficients. 
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Ratio of inaxinuin lift to miniuon 
profile drag. 
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Pig. 9 Increase of profile drag with lift coefficient. 



